Magnetic ordering of weakly coupled frustrated quantum spin chains 
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The ordering temperature of a quasi-one-dimensional system, consisting of weakly interacting 
quantum spin- 1/2 chains with antiferromagnetic spin- frustrating couplings (or zig-zag ladder) is cal- 
culated. The results show that a quantum critical point between two phases of the one-dimensional 
subsystem plays a crucial role. If the one-dimensional subsystem is in the antiferromagnetic-like 
phase in the ground state, similar to the phase of a spin chain without frustration, weak couplings 
yield magnetic ordering of the Neel type. For intra-chain spin-frustrating interactions larger than 
the critical one (at which the quantum phase transition takes place) , the quasi-one-dimensional spin 
system manifests a spiral magnetic incommensurate ordering. The obtained results of our quantum 
theory are compared with the quasi-classical approximations. The calculated features of magnetic 
ordering are expected to be generic for weakly coupled quantum spin chains with gapless excitations 
and spin-frustrating nearest and next-nearest neighbor interactions. 

PACS numbers: 75.10.Pq, 75.10.-b, 75.40.-s 



The interest in quasi-one-dimcnsional (quasi-lD) 
quantum spin systems has grown considerably during the 
last decades. The characteristic feature of quasi-lD mag- 
nets is the strong spin-spin interaction along one space 
direction, much stronger than the couplings along all 
other directions. The interest is motivated, on the one 
hand, by the progress in preparation of substances with 
well defined ID subsystems. Another reason for studying 
quasi- ID spin systems is the possibility to compare ex- 
perimental data with results of non-perturbative theories 
for ID models. Also, such systems often manifest quan- 
tum phase transitions that take place in the ground state, 
and which are governed by other parameters than the 
temperature, like external magnetic field, pressure, con- 
centration of impurities (internal pressure), etc. From 
the experimental viewpoint, quasi-lD spin-1/2 systems 
differ from other magnets due to special features in the 
behavior of their characteristics. For example, the tem- 
perature dependence of the magnetic susceptibility and 
the specific heat of quasi-lD spin systems with dominant 
nearest neighbor (NN) interactions reveal maxima (in a 
small external magnetic field at temperatures of the order 
of the exchange coupling constant along the distinguished 
direction) i For spin systems, which ID subsystems have 
a gapless spectrum of low-lying excitations at tempera- 
tures, much lower than the temperature of the maximum 
of the T-dependence of the susceptibility, the latter and 
the specific heat often manifest peculiarities, characteris- 
tic for phase transitions to low-temperature magnetically 
ordered states at critical temperatures. 

Also, a great attention has been recently given to spin 
systems with a spin frustration. For most of antiferro- 
magnetic (AF) systems the ground state corresponds to 
Neel-likc configurations. The standard quasi-classical de- 
scription of AF systems uses quantization of small devi- 



ations of vectors of magnetizations (magnetic order pa- 
rameters) of AF sublattices from their steady-state con- 
figuration. However, such a description of AF systems 
can be used for bipartite magnetic structures. For AF 
systems with a spin frustration the competing interac- 
tions produce a very high degeneracy of such steady-state 
configurations. Therefore, in most cases it is hopeless to 
use the approximation of magnetic sublattices. From the 
theoretical viewpoint the situation becomes even worse 
in quasi-lD spin systems with a spin frustration. For 
those systems quantum fluctuations are enhanced due to 
peculiarities in the ID density of states. This is why, 
according to the famous Mermin- Wagner theorem^ ID 
spin systems with short-range interactions cannot have 
any magnetic ordering even at T = 0. Thus, approxi- 
mate methods of theoretical physics often produce sig- 
nificant errors in the description of such systems. Hence, 
it is necessary to study them non-perturbatively, better 
exactly, which is, fortunately, possible for few cases for 
ID quantum spin systems. Probably one of the simplest 
and most known examples of a quantum spin system with 
a spin frustration is the Hciscnberg spin-1/2 chain with 
AF nearest neighbor (NN) and AF next-nearest-neighbor 
(NNN) interactions. The Hamiltonian of such a model 
can be written as 

(1) 

n n 

where J\ and J2 are the couplings between NN and 
NNN, respectively (here we consider only the case with 
an even number of spins N). Such a system is equiva- 
lent to a zig-zag spin ladder with obvious re-notation of 
indices. The system with the Hamiltonian ([1]) is, obvi- 
ously, spin-frustrated. Several limiting cases are known 
exactly. Namely, for J 2 = (or for J\ = 0) the Hamil- 
tonian |T|) is reduced to the Hamiltonian of one (or two 
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decoupled) Heiscnbcrg AF spin-1/2 chain(s). The ground 
state for those cases is a non-degenerate singlet, without 
long-range ordcrings, and the low-energy excitations are 
gapless spinonsji The other limiting case, for which the 
ground state is known exactly, is the so-called Majumdar- 
Ghosh point, J\ = 2J 2 .— In that case the Hamiltonian 
of the zig-zag spin ladder can be re- written as H.n nn = 
(^i/4) £„(S„ + S„ + i + S n+2 ) 2 -9A/4. The ground state 
is given by two degenerate singlets of the resonance va- 
lence bond type, without long-range ordering, and the 
low-lying excitations are gapped. For other values of 
the coupling constants it is, unfortunately, impossible, 
to obtain exact answers. Nevertheless, an approximate 
bosonization description and numerical calculations sug- 
gest that there is no long-range magnetic ordering in the 
system, and that for J 2 > 0.2411... J\ a spin gap is opened 
for the low-lying excitations.— A quasi-classical approxi- 
mation of the model yields the following. If one replaces 
the spin operators by classical vectors, two steady-state 
configurations are possible. The first one is the period 2 
commensurate and collinear AF Necl configuration, sta- 
ble for Ji > 4J 2 . The second one, stable for J\ < 4J 2 , is 
a noncollinear incommensurate spiral magnetic structure 
with the pitch angle cos cf> = — Ji/4J 2 . Such a description 
implies a long-range magnetic order. This means that it 
might be approximately valid for, e.g., a quasi-lD spin 
system, consisting of weakly coupled ID spin chains with 
NN and NNN AF interactions, at temperatures, lower 
than the ordering temperature. However, for quasi-lD 
spin systems with a spin gap A sp for low-energy excita- 
tions of their ID subsystems, a weak inter-chain coupling, 
as a rule, does not produce a magnetic ordering^ (This 
is plausible at least for spin systems with isotropic ex- 
change interactions: the exponential decay of the long- 
range spin-spin correlation function ex exp(— £/n) with 
a finite coherence length £ = hv/TrA sp , where v is the 
Fermi velocity of the low-lying spin excitations, conflicts 
with the magnetic order requiring asymptotically nondc- 
caying spin-spin correlation functions.) 

On the other hand, it is clear that spin frustration in a 
ID subsystem has to yield features in transitions to pos- 
sible magnetically ordered state for a quasi-lD system. 
Moreover, as follows from Ref. (see also Ref. 0) despite 
the fact that for most of studied compounds exchange 
constants satisfy the condition J 2 > 0.2411... Ji, the spin 
gap was not confirmed experimentally. To describe such 
experimental situations (i.e. quasi-lD AF spin systems 
with spin frustration of intra-chain interactions without 
a spin gap and with a weak inter-chain coupling), we 
consider another model, the Hamiltonian of which con- 
sists oiTtNNN with multi-spin interaction. Such a model 
is known to have gapless low-energy excitations. Those 
multi-spin interactions do not change the spin frustra- 
tion property from the classical viewpoint 4 The advan- 
tage of the proposed model is the exact integrability: The 
model permits an exact solution by means of the Bcthc's 
ansatz. We do not state here, naturally, that the model 
describes all features of the materials of current experi- 



mental interest^ However, many properties of the model 
are similar to what was observed in Ref. 0. At least, for 
this model the low-lying excitations are gapless. Hence, 
from this viewpoint, they qualitatively agree with the 
data of experiments, unlike the model with the Hamil- 
tonian H.NNN- Multiple spin exchange interactions are 
often present in oxides of transition metals, where a di- 
rect exchange between magnetic ions is complimented by 
a supcrexchange between magnetic ions via nonmagnetic 
ones^ Models with multi-spin interactions are believed to 
be closer to real quasi-lD magnets compared to standard 
ones with only NN spin couplings^ Multi-spin exchange 
models were introduced by Thouless already in 1965J^ 
Later similar models were used to study some cuprates^i 
and spin ladders For the consistent explanation of sev- 
eral experiments^ by means of inelastic neutron scat- 
tering, optical conductivity and nuclear magnetic reso- 
nance one needs to account for relatively large values of 
NNN spin-spin interactions and multi-spin interactions 
between four neighboring sites of the spin ladder (the 
so-called ring exchange). Similar four-spin interactions 
were used recently in the theory of 2D quantum spin sys- 
tems, where they regulate the quantum phase transition 
between the Neel-like ground state and the resonance va- 
lence bond solid oneJ^ The Hamiltonian of the ID sub- 
system of the quasi-lD model, studied in our work, has 
the form 

)(S 

n 

-(S„_iS n+2 )(S )) • (2) 

The model is also spin-frustrated. The classical coun- 
terpart of the model (if one replaces the spin operators 
by classical vectors) reveals a long-range magnetic or- 
dering with a Neel steady-state configuration, or with a 
spiral magnetic structure, where the four-spin ring ex- 
change renormalizes the spiral pitch angle as cos (f> = 
— 2Ji/(8J 2 — J4). However, quantum properties of the 
model with the Hamiltonian Hid differ from the one with 
TLnnn in a much more drastic way than of their classical 
counterparts. This can be seen from the exact solution 
(the exactly solvable model was introduced in Ref. [HI), 
which is known for the parametrization of coupling con- 
stants Ji = J(l — x), J 2 = Jx/2, J4 = 2Jx for any J 
and x (in what follows we shall consider J > 0, x > 0). 
For x = the model describes the Heisenberg spin-1/2 
AF chain. As one can see from exact results, the high de- 
generacy of low-energy states, caused by the spin frustra- 
tion of NN and NNN interactions, is removed by adding 
the ring exchange, which is also spin-frustrated. Accord- 
ing to the exactly known properties, the ground state 
of the model Eq. @ depends on the values of the cou- 
pling constant x and an external magnetic field H J£ For 
large values of H the model is in the spin-saturated (fer- 
romagnetic) phase. This phase has a trivial long-range 
magnetic order and gapped low- lying excitations. It is 
divided from other phases by the line of the second order 
quantum phase transition. For low values of x and H 
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FIG. 1: The phase diagram H — x of the one- dimensional 
integrable spin model. At the lines H s and H c second or- 
der quantum phase transitions to the ferromagnetic (spin- 
saturated) phase, and the ferrimagnetic spiral one, respec- 
tively, take place. In the point H = 0, x = x cr the second 
order quantum phase transition takes place. At the line H = 
for x > x cr the first order phase transition takes place. 



the model is in the phase, which behavior is similar in 
properties to the phase of the Heisenberg spin- 1/2 AF 
chain in a weak magnetic fieldi Notice that the model 
is in this phase for x < x cr = 4/tt 2 even at H = 0. At 
second order quantum phase transition takes place. 
For x > x cr and H = the model reveals an incom- 
mensurate ordering with nonzero spontaneous magneti- 
zation. For nonzero values of h and large enough values 
of x the model is in the incommensurate magnetic phase 
with nonzero weak magnetization^ The degeneracy in 
the direction of the spontaneous magnetization can be 
removed, if one first puts the system into an external 
magnetic field, and, then, removes the field. These two 
phases are divided from each other by the line of the 
second order quantum phase transition (quantum criti- 
cal point). Fig. 1 shows the ground state phase diagram 



of the one-dimensional integrable model. The mentioned 
quantum phase transitions can be observed in the tem- 
perature behavior of thermodynamic characteristics of 
the model, like the magnetic susceptibility and the spe- 
cific heat, that were also calculated exactly 

The goal of our present study is to find how the 
weak coupling between frustrated spin chains can pro- 
duce magnetic orderings, and what are specific features 
of such orderings. 

According to the above, one can suppose two differ- 
ent types of low-temperature magnetic ordering in the 
quasi-lD system under consideration. For the first one, 
the Neel ordering, one can write the magnetization of 
the n-th site, e.g., as M„ = Me, + (— l) n rriNe x , where 
&x,z are the unit vectors in the x- or z directions, M 
is the average magnetization, and tojv is the staggered 
magnetization in the direction, perpendicular to the ex- 
ternal field. Another possibility is the 3D generaliza- 
tion of the ground state incommensurate phase of the 
ID subsystem, the spiral incommensurate state with 
the magnetization of the n-th site, written as M„ = 
Me z + cos{Qn)m sp e z , where m sp is the modulated com- 
ponent of the ^-projection of the magnetization around 
the average magnetization M, and Q = tt(1 — 2M) is the 
wave vector of the ID modulated structure. 

As usually, we study the weak inter-chain coupling J' 
in the mean field approximation. In that approximation 
in the Neel phase we can write the mean field Hamilto- 
nian of the ID subsystem as 

H™ f = -Hid -{H- zJ'M) ^ S z n 

n 

-h N £](-l) n S£ + const , (3) 

n 

where h n = zJ'ttin, and z is the coordination number. 
For the spiral phase the mean field Hamiltonian is 

H™J = Hid - {H - zJ'M) ^ S z n 

n 

—h sp 2^ cos(Qn)S z l + const , (4) 

where h sp = zJ'm sp . Renormalization group-like 
approach^ implies that both and h sp are relevant per- 
turbations. They generate spin gaps AEjy ~ /i^ 4 ^ 

and AE sp ~ /i^p ^ , respectively, for low-energy ex- 
citations. Here rj is the correlation function exponent, 
see below, which determines the asymptotical behavior 
of the spin-spin correlation functions of the ID subsys- 
tem in the conformal limitiii 

The order parameters m^r and m sp (or and h sp ) 
have to be determined self-consistently. In the mean field 
approximation the corresponding self-consistency equa- 
tions can be written as 

m N ,s P = M NtSp (H, h N:Spi T) , (5) 

where MN tSp (H, h^^ sp ,T) is the magnetization per site 
of the ID subsystem in the effective field H — MzJ' and 
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h,N,sp a t the temperature T. Then the transition tem- 
perature to the ordered state has to be determined from 
the equation 



1 = zJ'xN,st , 

'dM NtSt (h,h N , sp ,T) 



XN.st 



dh 



N,sp 



(6) 



The non-uniform susceptibilities of the ID subsystem at 
low temperature can be found as 

Xa (q,T) = -iY / I dte-^ n e(t)([S a (n,t),S a (0 7 0)}) T , 

(7) 

where q is the wave vector, a = x, z, and (...)t denotes 
the thermal average at the temperature T. Asymptotics 
of the correlation functions for an integrable spin chain 
can be obtained in the conformal field theory limit X For 
the model with the Hamiltonian Tim it was done in 
Ref. H?l and we can write for the ground state correlation 
functions 



M 



(s x n s x ) « (-1)' 



B* cos(Qn) 
[n 2 - {vt) 2 } e * 
C 

[n 2 - (vt) 2 ] 9 ± 



(8) 



where v is the Fermi velocity of low-energy excitations, 
6 Z = Z 2 {= 1/2t/), 9± = l/AZ 2 = (77/2), Z is the dressed 
charge of low-lying excitations, B* , and C are nonuni- 
versal constants. In particular, we see, that the sym- 
metry of the ground state is lower than the symmetry 
of the Hamiltonian, caused by the ordering, i.e. for our 
model one deals with the manifestation of the Goldstonc 
theorem. Eqs. |jHJ) can be extended for weak nonzero 
temperatures using the conformal mapping (n — vt) — > 
(2v/ttT) sinh[7rT(n— vt)]X Then, calculating susceptibili- 
ties according to Eq. |7]) (we use the main approximation) 
for q = ir for the Neel phase and for q = Q for the spiral 
incommensurate phase, we obtain the expressions for the 
ordering temperatures 



T, 



N 



V 

2^ 



<4-n(H)W2, 



n 



(9) 



and 



T — 



V 

2~ 



B* 



zJ' 



— 1 B 2 ( — 

2 V J \4q' 



2<q - 1 



v \2r) J V 4r 7 ^77 

(10) 

where B(x, y) is the Euler's beta function. In those ex- 
pressions the Fermi velocity and the critical exponent 
77 (or the dressed charge Z) can be calculated using the 
Bethe ansatz results^ Then, the question to be answered 
is: Which ordering temperature, Tjv or T sp , is higher for 
the quasi-lD spin chain with the spin frustration? 

In what follows we limit ourselves with the case H = 
for simplicity. In this situation the effective Fermi veloc- 
ity can be written as v — (n/2)J\l — (x/x cr )\M- Consider 



first the ground state phase x < x cr for a ID subsystem, 
which is similar to the ground state of the Heisenberg 
spin-1/2 chain with only NN AF interactions. In this 
case we have M = and ?y = lJ2 For this case we can 
use B* = C ~ 0.2fi£ One can see that in this case (i.e. 
rj = 1) T N = T sp = {CzJ'/2tt)B 2 (1/4:, 1/2). We see 
that the critical temperature does not depend on J and 
x (obviously, any nonzero magnetic field h ^ 0, or an in- 
clusion of a magnetic anisotropy will change this result). 
To get the J- and x-dependences even for H = and for 
the magnetically isotropic case, one has to include log- 
arithmic corrections^ reproducing the known result for 
x = 0l2. 



^v = ^B 2 (l/4 J l/2) 1 



7T 2 J|1 - (x/x cr )\ 



'hi 



CzJ'B 2 (l/A, 1/2) 



(11) 

which is valid, naturally, when the argument of the log- 
arithm is larger or equal to 1. Next, let us consider 
the ground state of the ID subsystem for x > x cr , 
which ground state has a spontaneous magnetic order- 
ing. This spontaneous magnetization M(x > x cr ) 7^ 
is connected with holes in the ground state distribution 
of quantum numbers, called rapidities, which form the 
Dirac sea of the ID subsystem! 15 i 17 Those holes appear 
only for x > x cr i 15 i 16 i 17 Notice that in the previous case, 



there arc no holes in the Dirac sea, and the 



ground state rapidities can have any value in the range 
It is impossible to find an analytic 
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solution for 77 in this case. We see that rj = 1 + a, with 
< a < a max < 1 when 1 < x/x cr < 00. Unfortu- 
nately, in this case we cannot obtain the values of the 
non-universal constants. We can only suppose that they 
are also of the order of 0.05 — 0.2J^ It is easy to see that 
for most of the values of J, J' and x the temperature 
of the transition to the spiral incommensurate state for 
x > x cr is higher than the Neel temperature. In Fig. 2 
we plotted Tjv (lower surface) and T sp (upper surface), 
for J = 1, J' = 0.01, z = 4 and B = C = 0.2 as functions 
of a and y = (x/x cr ). It turns out that a is a function 
of x/x cr also, but, unfortunately, one cannot find this 
dependence analytically. Notice that for a = (77 = 1, 
x = x cr ) both critical temperatures coincide. The Neel 
temperature can be larger than the temperature of the 
transition to the spiral incommensurate phase only in the 
vicinity of the quantum phase transition x = x cr for 77 
being very large (close to 2, which seems to be an over- 
estimation, cf. ReflTBT). For all other values of x > x cr we 
get T sp > TV- Hence, we can conclude that for a quasi- 
lD system, consisting of weakly coupled spin-1/2 chains 
with AF spin-frustrating NN and NNN interactions and 
with the four-spin ring exchange, the low-temperature 
ordering depends on the behavior of ID subsystems. For 
small values of the NNN interaction, the quasi-lD sys- 
tem undergoes a transition to the magnetically ordered 
AF Neel state. On the other hand, if the exchange con- 
stant of the NNN interactions exceeds the critical value, 
at which a quantum phase transition to the incommcn- 
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FIG. 2: The Neel temperature (the lower surface) and the 
ordering temperature of the transition to the spiral incom- 
mensurate phase for a quasi-lD spin-1/2 chain with the spin 
frustration, caused by AF NN, NNN interactions and the ring 
exchange as functions of the deviation of the critical exponent 
a — 1 — 77 and y — x/x cr , which shows how close the quantum 
critical point y = 1 is. 

suratc state with the weak spontaneous magnetization 
takes place, a weak coupling between ID subsystems pro- 
duces the transition to the magnetically ordered incom- 
mensurate spiral state. The ordering temperature in the 
latter is 

J|l-2/i r D . 2zJ' , / 7T \ 

i,„ = B : sin — x 

P 4 [ 7rJ|l-y| \2{l + a)) 

B2 {w^rw^))} ■ (12) 

Our quantum analysis qualitatively agrees with the 
quasi-classical description of the considered system. 
However, a difference between the quantum system and 
its quasi-classical counterpart exists: In the quantum sys- 
tem magnetic ordering takes place for J' 7^ 0, only. We 
expect analogous expressions to Eqs. (11,12) to be valid 
also for the other critical point between a ferromagnetic 
and a spiral phase, i.e. if J\ < 0. In this context, the 
determination of the magnetic order below their low- 
temperature phase transitions at few K for Li 2 ZrCu04 
and Pb[Cu(SO) 4 (OH) 2 ],— both being close to that crit- 
ical point, would be of interest. 

We expect that an inter-chain ferromagnetic 
interaction,— for x < x cr has to produce the AF 
low-temperature ordering of a special type. Namely, we 
expect low temperature ordered phase, which consists 
of ferromagnetic planes with alternating magnetizations 
of planes. On the other hand, for x > x cri a ferro- 



magnetic inter-chain interaction is expected to produce 
ferrimagnetic low-temperature ordering with a nonzero 
spontaneous total magnetization. 

Finally, let us consider what happens, if one studies 
the situation with only NN and NNN couplings, with- 
out the ring exchange. In that case, for J 2 < 0.2411...Ji 
the quasi-lD system undergoes a phase transition to the 
Neel state, due to weak couplings between chains. For 
J 2 > 0.2411...Ji, a spin gap is opened for low- lying ex- 
citations of the ID subsystem, and the ordering tem- 
perature goes to zero. This case seems to contradict 
known experiments, in which magnetic ordering was ob- 
served even for J 2 > 0.2411... Ji^ Therefore, we can con- 
clude, that for some real compounds with the properties 
of quasi-lD spin systems with spin-frustrating interac- 
tions in their ID subsystems some additional spin-spin 
interactions, like the ring exchange, studied in this pa- 
per, probably exist, which close the spin gap and give 
rise to magnetic orderings at low temperature. 

In summary, the ordering temperature of a quasi- 
one-dimensional system, consisting of weakly interacting 
quantum spin-1/2 chains with antifcrromagnctic spin- 
frustrating couplings (or zig-zag spin ladder) is calcu- 
lated. Our results show that the quantum critical point 
between the two phases of the ID-subsystem plays an im- 
portant role. If the one-dimensional subsystem is in the 
ground state in an antiferromagnetic-like phase, similar 
to the phase of a spin chain without frustration, weak 
couplings yield a magnetic ordering of the Neel type. On 
the other hand, for intra-chain spin-frustrating interac- 
tions larger than the critical one (at which the quantum 
phase transition takes place), an incommensurate spi- 
ral magnetic ordering of the quasi-one-dimensional spin 
system takes place. The obtained results of the quan- 
tum theory are compared with the quasi-classical ap- 
proximations. We expect that the calculated features 
of the magnetic ordering are generic for weakly coupled 
quantum spin chains with gapless excitations and with 
spin-frustrating nearest and next-nearest neighbor inter- 
actions. While up to now we do not know quasi-one- 
dimensional systems with NN AF spin interactions and 
large AF NNN ones in the spiral phase at low temper- 
atures (cf. Ref. H, see, though, Refs. I22II23I ), we believe 
that our results can be used for comparison with the ob- 
served temperatures of magnetic orderings in other spin 
frustrated quasi-lD quantum spin systems. 
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